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Abstract.  In  recent  paper  [K.  Ito,  et  al.,  Jpn.  J.  Appl.  Phys.  40,  2558  (2001)],  we  present  a fast 
and  sufficiently  accurate  procedure  to  construct  the  potential  and  the  electric  field  distribution 
from  the  observed  electron  density  distribution  in  two-dimensional  plane.  Using  this  procedure,  we 
obtain  quantitative  agreement  between  (he  observed  velocity  of  the  clump  motion  and  the  calculated 
velocity  of  that  determined  by  the  field  analysis.  In  addition,  the  strong  deformation  of  the  clump  is 
observed  when  that  moves  in  the  shear  of  the  background  vortex. 


INTRODUCTION 

The  dynamics  of  a nonneutral  electron  plasma  that  is  strongly  magnetized  and  ho- 
mogeneous along  the  magnetic  field  is  equivalent  with  the  two-dimensional  (2D)  Euler 
fluid  (1,  2].  The  correspondence  between  the  electron  density  n(x,v)  and  the  vorticity 
^(.v,y)  and  that  between  the  electric  potential  q>(jc,  >>)  and  the  stream  function  v) 
eases  the  diagnostics  and  analyses  of  the  electron  vortex  motion  compared  with  those 
for  neutral  fluid  vortices.  The  2D  structures  of  the  density  (vorticity)  distribution  is  de- 
termined experimentally  in  terms  of  the  luminosity  distribution  of  electrons  collected  on 
a phosphor  screen  that  serves  also  as  a collector  of  the  total  charge  to  determine  the  total 
number  of  electrons  [3,  4]. 

From  the  observed  2D  density  distribution,  we  determine  the  electrostatic  potential 
and  electric  field  distribution  in  the  2D  plane  by  employing  Fourier-Bessel  series  expan- 
sion of  the  observed  density  distribution  for  fast  numerical  processing  [4],  The  potential 
distribution  $(;<:, y)  conforms  to  the  stream  function  vj/(-v,y)  in  terms  of  fluid  dynamics, 
so  that  many  aspect  of  fluid  dynamics  can  be  presented  in  the  language  of  electromag- 
netism. This  fast  numerical  analysis  is  essential  for  detailed  and  extensive  studies  of 
vortex  dynamics  or  turbulence. 

We  have  experimentally  observed  radial  oscillations  in  the  trajectory  of  a clump  of  a 
vorticity  moving  up  the  slope  of  the  background  vorticity  distribution  (3].  This  oscilla- 
tion is  not  expected  from  the  linear  theory  [5],  and  we  were  led  to  the  interpretation  that 
the  main  mechanism  responsible  for  the  oscillations  is  the  azimuthal  electric  field  associ- 
ated with  the  differential  rotation  of  the  perturbed  density  distribution  in  the  background 
behind  the  clump  [3]. 
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Using  this  field  analysis,  we  have  obtained  quantitative  agreement  between  the  ob- 
served speed  of  the  oscillation  and  the  calculated  speed  determined  by  the  field  analysis. 
It  also  demonstrates  the  influence  of  structures  in  the  background  vortex  on  the  dynam- 
ics of  coherent  vortices  [6,  7,  8,  9].  On  the  other  hand,  it  is  observed  that  the  clump  is 
stretched  in  the  background  vortex  shear  along  the  stream  function.  There  is  a possibility 
that  it  also  influences  the  vortex  motion. 

The  experiment  is  carried  out  with  the  Malmberg  trap.  The  details  of  the  experimental 
configuration  have  been  reported  in  refs.  3 and  4. 


DERIVATION  OF  POTENTIAL  FIELD  FROM  OBSERVED 
DENSITY  DISTRIBUTION^ 


In  the  Malmberg  trap  configuration,  the  magnetic  field  B o is  homogeneous  and  the 
axial  distribution  of  the  vacuum-field  potential  <j> v is  fiat  in  the  trap  region.  Therefore,  we 
can  safely  approximate  that  the  chord-averaged  density  n(x,y)  — n i{x,y)/L  represents 
the  actual  density  distribution  of  the  electrons.  Here,  L is  the  axial  length  of  the  trap 
region  expected  from  <j>v.  Because  L is  much  larger  than  the  transverse  dimension  of 
the  electron  density  distribution,  we  can  introduce  a 2D  analysis  for  determining  the 
potential  distribution  in  the  main  part  of  the  trapped  region  by  neglecting  the  effects 
from  the  ends  [4]. 

The  potential  distribution  generated  by  the  electrons  is  given  in  a 2D  form  <()(r,  0)  that 
satisfies  the  Poisson’s  equation. 


fld_  £ 

\ r dr  dr  + r2  920  ) 


<Kr>0)  = ~~n(r,Q). 
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The  contribution  of  the  conducting  wall  surrounding  the  plasma  enters  as  a boundary 
condition  <(>(/?,  0)  = 0 at  the  wall  of  a radius  R. 

First,  we  Fourier-expand  Poisson’s  equation  in  the  azimuthal  angle  0 to  obtain  a set  of 
ordinary  differential  equations  in  the  radial  coordinate  r.  <b  f/(r)  and  nCm(r)  are  the  wth 
Fourier  component  of  the  potential  and  the  density,  respectively.  The  density  part  n c„]s(r ) 
are  calculated  from  the  observed  CCD  images.  We  further  expand  n nf(r)  in  a series  of 
Bessel  functions  Jm(z)  to  determine  the  associated  terms  of  the  potential  that  is  obtained 
as, 


*M>  = 
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where  Jm(Xnm)  = 0 with  n = 1,2,3,....,  and 
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Here  the  upper  limit  of  the  azimuthal  mode  number  mmax  and  the  radial  mode  num- 
ber nmax  are  taken  as  mmax  — 100  and  nmax  — 100.  The  calculation  of  the  potential  re- 
quires about  1 min.  on  a personal  computer.  Once  the  potential  distribution  <f>(r,  0)  is 
constructed,  the  electric  field  can  be  determined  from  <t>(r,  0). 

Similarly,  we  can  also  reconstruct  the  density  profile  if  we  remove  the  constant  ejn 
and  replace  the  coefficients  p c„f„  with  acmsn  — —(%,„„//?) 2 in  eq.  (2).  The  discrepancy 
Ae  between  the  original  density  profile  n„  and  the  reconstructed  profile  iifu  may  be 
defined  as 

A,  = (3) 

I n„ 

by  summing  all  the  pixels.  In  most  cases,  A<,  is  found  to  be  less  than  1%. 


DYNAMICS  OF  A CLUMP  IN  THE  BACKGROUND  VORTEX 

In  this  section,  we  show  an  experimentally  observed  dynamics  of  a clump  in  a 
background  vorticity  distribution.  In  our  experiments,  generally,  the  radius  of  the  clump 
is  very  small  compared  with  that  of  the  background  vortex,  and  the  circulation  T f of 
clump  is  less  than  a few  percent  of  T The  circulation  f is  proportional  to  the  electron 
number  N.  However,  the  vorticity  of  the  clump  C,c  is  ten  times  as  high  as  that  of  the 
background  vortex  Accordingly,  we  can  assume  the  clump  as  the  strong  point  vortex 
in  our  experiments. 

Figures  1 show  typical  time  evolution  of  the  2D  vorticity  distribution.  Here,  the 
electron  number  of  the  clump  and  background  vortex  is  Nc  = 1.2  x 1()7  and  Nb  = 
1.9  x 108,  respectively.  The  upper  panels  are  total  vorticity  distribution  observed 


FIGURE  1.  Time  evolution  of  2D  vorticity  distribution.  The  upper  panels  show  the  total  vorticity 
distribution  observed  with  CCD  camera  . The  lower  panels  show  the  perturbed  part  of  the  background 
vortex  from  the  initial  background  vortex  distribution. 
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with  CCD  camera.  It  is  observed  that  the  clump  climbs  up  the  background  vorticity 
gradient  as  it  rotates  creating  a spiral  streak  behind  it. 

The  lower  panels  show  the  perturbed  part  of  the  background  vortex  distribution 
8 C,b(x,y,t)  which  is  determined  by  subtracting  the  initial  background  distribution 
^(jc,y,0)  from  each  distribution  C,(x,y,t)  in  the  upper  panels.  It  is  observed  that  the 
negative  (positive)  perturbation  in  front  (behind)  of  the  clump.  It  is  considered  that  this 
deformation  of  the  background  vortex  generates  the  radial  electric  field,  and  the  E x B 
drift  caused  by  this  field  drive  the  radial  motion  of  the  clump. 


ANALYSIS  OF  THE  CLUMP  DYNAMICS 

Figure  2(a)  shows  the  radial  distance  rc  of  a clump  from  the  center  of  the  background 
vorticity  against  the  time  with  different  circulation  of  clumps  (N  c=  (2.4  — 6.3)  x 106) 
when  the  background  vorticity  profile  is  fixed  (Nb  — 1.9  x 108).  The  clump  begins  its 
free  motion  at  15  rs  on  disconnection  from  the  electron  source.  Solid  lines  represents 
radial  distance  calculated  from  the  Schecter’s  theoretical  model  [5]  using  the  full  ex- 
perimental parameters  of  the  initial  distribution.  The  experimentally  observed  orbit  of 
the  clump  is  in  almost  agreement  with  the  orbit  obtained  from  the  theoretical  model.  In 
the  experiment,  however,  we  can  observe  the  oscillation  of  the  orbit  of  the  clump  mo- 
tion. We  consider  that  the  oscillation  of  the  orbit  is  driven  by  the  azimuthal  electric  field 
caused  by  the  strong  deformation  of  the  background  vortex  [3]. 

Figure  2(b)  shows  the  radial  velocity  of  the  clump  motion  in  the  case  indicated  by 
open  triangles  in  Fig.  2(a).  The  open  triangles  represent  the  radial  velocity  calculated 
from  the  observed  radial  distance  of  the  clump.  The  closed  triangles  indicate  the  radial 
velocity  driven  by  E x B drift  when  the  electric  filed  acting  on  the  clump  is  calculated 
with  the  field  analysis.  The  solid  line  represents  the  radial  velocity  calculated  form  the 


FIGURE  2.  (a)  The  symbols  are  indicate  the  observed  radial  distance  rc  of  the  clump  from  the  center 

of  the  background  vortex  for  different  1VC/106  =2.4(0),  4, 1(A),  6. 3(D).  Solid  lines  represents  the  radial 
distance  calculated  from  the  Schecter’s  theoretical  model,  (b)  Radial  velocity  of  the  point  vortex  in  the 
case  Of/Vc=4.1xl06  in  Fig.  2(a).  Open  triangles  represent  the  directly  observed  velocity  from  observed 
vorticity  distribution  and  closed  triangles  indicate  the  velocity  calculated  with  the  field  analysis  method. 
Solid  line  shows  the  velocity  expected  from  the  theoretical  model. 
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FIGURE  3.  (a)  Typical  image  of  a clump  stretched  with  the  shear  of  the  background  vortex  at  t = 
30  fts.  Thin  lines  represent  the  stream  function.  Thick  line  indicates  the  1/e  contour  of  the  vorticity.  (b) 
Open  triangles  represent  the  ratio  of  the  major  axis  to  minor  axis.  The  solid  line  indicates  the  shear  of  the 
background  vortex. 


theoretical  model. 

The  radial  velocity  calculated  with  the  field  analysis  is  in  good  agreement  with  the 
observed  velocity.  It  is  noticed  that  the  reverse  of  the  radial  velocity  is  reproduced  from 
the  field  analysis.  Therefore,  we  actually  confirm  the  reverse  of  the  radial  velocity  driven 
by  the  E x B drift  by  the  field  analysis.  It  indicates  that  this  fast  numerical  analysis  is 
essentia]  for  detailed  and  extensive  studies  of  vortex  dynamics. 

It  is  observed  that  the  maximum  velocities  calculated  from  the  experimental  data  are 
slow  compared  with  the  theoretical  value.  We  also  consider  that  it  originates  in  a strong 
deformation  of  the  background  vortex  and  that  of  the  clump  mentioned  below. 

We  examine  the  influence  which  the  background  vortex  exerts  on  the  clump.  Figure 
3(a)  shows  the  typical  vorticity  distribution  around  the  clump  in  the  background  vortex. 
Thin  lines  represent  the  stream  function  and  thick  line  indicates  the  1/e  contour  of  the 
clump  vorticity.  It  is  observed  that  the  clump  is  stretched  along  the  stream  function. 

Figure  3(b)  plots  the  ratio  a of  the  major  axis  to  minor  axis  of  the  clump  (A)  and  the 
rotation  shear  A = - rdQ.b/dr  of  the  background  vortex  (solid  line)  as  a functions  of 
the  time.  Here,  Q(,(r)  is  the  angular  frequency  of  rotation  of  the  background  vortex.  The 
clump  is  round  (aw  1),  immediately  after  the  free  motion  of  clump  starts.  After  that, 
the  clump  is  stretched  swiftly  along  the  stream  function.  The  clump  is  stretched  most 
(a  « 2.2)  at  / = 30  jts  when  the  shear  A of  the  background  vortex  is  strongest  and  the 
radial  velocity  of  the  clump  is  fastest  as  shown  in  Figs.  3(b)  and  2(b).  When  the  clump 
approaches  the  center  of  the  background  vortex,  the  shear  A of  the  background  vortex 
becomes  weak  and  the  clump  rounds  again.  It  is  confirmed  that  a clump  is  stretched  with 
a shear  of  a background  vortex. 

In  a numerical  simulation  of  neutral  fluid,  it  has  been  considered  that  the  deformation 
of  the  vortex  is  very  important  to  the  merging  between  vortices  [10].  We  have  reported 
that  the  merging  between  the  clumps  is  accelerated  in  the  background  vortex.  On  the 
other  hand,  we  have  observed  that  two  clumps  continue  to  revolve  around  each  other 
with  a fixed  distance  though  they  approach  in  close  distance  in  the  background  vor- 
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tex  [7].  We  consider  that  the  clump  dynamics  and  the  interaction  between  these  in  the 
background  vortex  are  controlled  by  not  only  the  transformation  of  the  background  vor- 
tex but  also  the  clump  deformation  caused  by  the  interaction  between  the  clump  and  the 
background  vortex. 


CONCLUSIONS 

We  have  developed  a fast  procedure  to  numerically  construct  the  2D  field  distribution 
from  the  observed  image  data  directly  related  to  the  electron  density  distribution.  We 
have  examined  the  dynamics  of  the  clump  in  the  background  vorticity  distribution  using 
this  field  analysis.  As  a result,  we  have  obtained  the  quantitative  agreement  between 
the  observed  velocity  of  the  clump  motion  and  the  calculated  velocity  determined  by 
the  field  analysis.  It  includes  the  oscillation  of  the  clump  orbit.  It  is  consider  that  this 
method  is  useful  for  extensive  studies  of  vortex  dynamics  or  turbulence. 

We  also  examine  the  deformation  of  the  clump.  It  is  observed  that  the  clump  which 
has  the  strong  vorticity  is  stretched  in  the  background  shear.  We  have  considered  that  the 
clump  deformation  caused  by  the  background  vortex  shear  plays  a role  as  important  as 
the  background  vortex  deformation  in  the  dynamics  of  a clump  and  interaction  between 
these  in  the  background  vortex. 
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